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We compare the charge transport characteristics of heavy doped p ++ - and n ++ -Si-alkyl chain/Hg 
junctions. Photoelectron spectroscopy (UPS, IPES and XPS) results for the molecule-Si band 
alignment at equilibrium show the Fermi level to LUMO energy difference to be much smaller than 
the corresponding Fermi level to HOMO one. This result supports the conclusion we reach, based on 
negative differential resistance in an analogous semiconductor-inorganic insulator/metal junction, 
that for both p ++ - and n ++ -type junctions the energy difference between the Fermi level and 
LUMO, i.e., electron tunneling, controls charge transport. The Fermi level-LUMO energy difference, 
experimentally determined by IPES, agrees with the non-resonant tunneling barrier height deduced 
from the exponential length-attenuation of the current. 



I. INTRODUCTION 

Molecular electronics describes charge transport pro- 
cesses whereby molecules serve as active elements (e.g., 
rectifiers, switches, sensors) or passive ones (resistors or 
surface passivating agents) in electronic devices^. Since 
the emergence of this field it has expanded greatly and 
now includes several types of device configuration (e.g., 
two-terminal junctions^', three-terminal junctions^ and 
electrochemical devices^), where the molecules serve as 
direct current carriers ^ or in directly affect the electrical 
properties of a junctiori 11 * 12 !. 

Despite the diversity of investigated systems, funda- 
mental questions regarding the mechanisms for electri- 
cal current passing through molecules between two elec- 
trodes and the possibility of gaining predictive power and 
control over the electrical properties of molecular junc- 
tions remain mainly unsolved. Detailed discussions on 
the limitations of existing transport models can be found 
elsewhereP^HHI i n short, even for one of the simplest sys- 
tems, i.e. that of alkyl chains between Au electrodes, 
there is a large discrepancy between the average tunnel 
barrier height, extracted f rom c urrent - voltage (7 — V) 
measurements (~ 1.2eV)P^^ and the barrier that 
is expected from the experimentally determined elec- 
trode work function and alkyl monolayer ionization po- 
tential, i.e. the energy difference between the Fermi 
level and highest occupied molecular orbital (HOMO) as 
found by Ultraviolet Photoelectron Spectroscopy, UPS 
(~ 5 eV)P^In spite of this apparent discrepancy, HOMO- 
dominated transpor t is the prevailing concept for such 
j unctions EH t 16 l 23 EH l n the present study we challenge this 



concept and find that transport is dominated by the en- 
ergy difference between the Fermi level and the lowest 
unoccupied molecular orbital (LUMO) of the molecules. 

It would seem that having a semiconductor (SC) in- 
stead of a metal as one of the electrodes in a molecu- 
lar junction increases the difficulty in understanding the 
mechanisms that govern charge transport through the 
junction. Indeed, with moderately doped SCs, the SC 
depletion layer has a large effect on the overall observed 
current and if the band bending is very large, it com- 
pletely overwhelms any molecular effect.^ As we shall 
see, though, the intrinsic asymmetry of a junction with a 
SC, with respect to charge carrier type, provides unique 
information on the energy levels involved in transport, 
in ways that are impossible when using only metal elec- 
trodes. 

To find how charge is transported through a mono- 
layer of saturated alkyl chains we measure and analyze 
electrical transport across heavy-doped semiconductor - 
molecular insulator / metal junctions. Our results pro- 
vide unequivocal evidence that charge transport in these 
junctions is controlled by the energy difference between 
the metal Fermi level and the LUMO of the molecule, 
\Ep — LUMO |, rather than by the energy difference with 
the HOMO, \E F - HOMO|. This finding agrees with 
complementary spectroscopic measurements of the sam- 
ples without Hg contact, but challenges the pre vailing 
concept of HOMO-dominated transport jHl t 16 l 23 l 24 l Highly 
doped n ++ - and p ++ - Si-C„H2 n +i/Hg junctions (n — 
2, 14, 16, 18) were prepared by alkylation of freshly etched 
P- or B-doped (AT d/a w 10 19 c m - 3 ) Si(l 1 1) surfaces. 
Detailed descriptions of the preparation and character- 
ization of both 'long' (C14-C18) and 'short' alkyl chain 



1.0- 


a 


i i i i i i i i i i 
b 


l . < . . l < . . < 1 i i 
C / 






\ 


/' 
// 


0.8- 




\ 


// 












V 


ri 
1 / 


0.6- 




\ 


it 
II 






\ 


jl 


0.4- 


\pH-+ v m-+ 


Jt 
ft 

ft 






\ 


fi 






\ 


Ji 


0.2- 




\ 


fi 
// 






\ 




0.0- 




; 





-1.5 -1.0 -0.5 0.0 0.5 1.0 



Metal bias (V) 

FIG. 1: Normalized G-V curves for p ++ - (solid line) and n ++ - 
(dashed line) Si-Ci6H33/Hg junctions. The conductance of 
the n ++ -junction is parabolic around V, while for the p ++ - 
junction at negative bias on the metal there is a voltage range, 
where the conductance increases slowly (shaded area), fol- 
lowed by a sharp increase. The transition bias between the 
above-mentioned regions is ~ 1.1 V, very close to the Si band 
gap (~ 1.12 eV). (a), (b) & (c) refer to voltage ranges, shown 
in Fig. 2 and discussed in the text, for the p ++ junction only. 



(C2) samples can be found elsewhere, and ar e given in 
the Suppl. Inf. (section 1 and section 2)! 2S | 27 1 



II. RESULTS AND DISCUSSION 

In Fig. n] the normalized conductance of p ++ -/n ++ - 
Si-Ci6H33/Hg junctions is compared on a linear scale. 
While conductance for the n ++ junction (dash) is 
parabolic around V (i.e. similar to the behaviour 
expected for Metal -Insulator -Metal (MIM) tunnel 
junctions^, for the p ++ junction (solid) there is a neg- 
ative bias range where conductance increases slowly 
(shaded), followed by a sharp increase at ~ 1.1 V, which 
closely corresponds to the Si band gap (1.12 eVp*H 

This behaviour is highly reproducible and indepen- 
dent of the molecular length (C14H29-C18H37; see ref. 
|2"T1 for J — V curves of n ++ -Si data and Fig. [5] below 
for the p ++ -Si ones). The results of the p ++ junction 
at negative applied bias can be understood by consid- 
ering work on Metal-inorganic Insulator- (near) degen- 
erate Semiconductor (MIS) tunnel diodes by Esaki, Sze, 
Dahlke and co workers J22H22] which forms the basis for 
the present study. Schematic band diagrams of the rel- 
evant model are shown in Fig. [2] (top), along with an 
experimentally measured semi log J — V curve of a typ- 
ical p ++ Si-Ci6H 3 3/Hg junction. We identify the undu- 
lating behaviour at low to moderate negative bias (re- 
gion b of Fig. [2j which is equivalent to the shaded region 
111 Fig. [TJ with a region where negative differential resis- 
tance (NDR) was predicted originally for the MIS system 
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Metal bias (V) 

FIG. 2: Schematic band energy diagrams for MIS structures 
with p ++ semiconductor substrates and LUMO-dominated 
transport (Top) and semi- log J-V curve of a typical p ++ Si- 
Ci6H33/Hg junction (averaged over 12 scans) (Bottom) with 
the (a), (b) and (c) labels corresponding to the diagrams 
above. Three major bias regimes are considered in the dif- 
ferent energy schemes (top) and separated in the J — V 
plot by vertical dotted lines (bottom), corresponding to a) 
metal Fermi above (i.e., closer to vacuum level) the SC 
gap (V < —(Eg + C))i b) metal Fermi within the SC gap 
(-(E s + C) < V < — C); and c) metal Fermi below the SC 
gap (V > —( ), where E s and £ are the forbidden gap and 
the energy difference between the SC Fermi level and the va- 
lence band edge. TOP: tj>t is the energy difference between the 
tunneling electrons (marked by a horizontal arrow) and the 
molecular level (spatially averaged), which is taken as the es- 
timated effective tunneling barrier height. The insulator level 
is assumed to vary across the insulator. Electron tunneling 
into the VB (as in b) is practically identical to holes tunneling 
in the opposite direction. Surface states, band-bending and 
image forces are neglected for simplicity. 



studied by Esaki et al. That effect is smeared out in the 
J — V curve, because of the effect of interface states, an 
effect that, together with band bending and image forces, 
is, for simplicity's sake, neglected in the schematic band 
diagrams. In the figure, only E-p -LUMO (i.e., electron 
tunneling) dominated transport is considered (LUMO is 
depicted by thick, black line) and other cases are dis- 
cussed below. 

As long as the metal Fermi level does not align with 
the SC forbidden energy gap (i.e., regions a, c) transport 
is completely equivalent to tunneling in the more com- 
mon Metal-Insulator-Metal (MIM) tunnel junctionPEH 
In regions a and c, the bias acts to reduce the averaged 
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tunneling barrier by ~ |V|/2. With a moderate nega- 
tive voltage to the metal (Fig. [2| the metal Fermi level 
moves to energies that are in the forbidden band gap 
of the SC. For such a case, neglecting interface states- 
assisted tunneling (IS AT), Esaki and Stiles predicted an 
effective increase in the tunnel barrier height 4>t (i-e., 4>t 
(b) > <f> t (a,c)), and NDR, as has indeed been observed 
experimentally. 32 If ISAT is not neglected, instead of 
NDR a ln(J) — V plot will show undulating behaviour 
(as we show in Fig. [2] voltage range (b), for molecular 
junctions), which is strongly influenced by the interface 
state density and distribution in the band gapP^ In our 
case ISAT has a clear influence on the J — V data and, 
therefore, we observe undulating J— V behaviour, instead 
of pure NDR. Undulating J — V occurs because in region 
b the highest energy electrons have no matching states to 
tunnel into, and transport proceeds via deeper electrons, 
aligned with the SC VB. For those electrons the barrier 
is now increasing as ~ |V]/2 instead of decreasing with 
|V|, as in band-to-band tunneling (ranges a,c). In con- 
trast to p ++ junctions, n ++ junctions are not predicted 
to show undulating J — V plots for LUMO-dominated 
tunneling because the SC Fermi and the tunneling levels 
are on the same side of the forbidden SC gap, regardless 
of applied bias (see Fig. S3. lb of Suppl. Inf.). Simply 
put, undulating J — V is predicted to occur if tunneling 
proceeds via carriers of opposite type to the SC majority 
carriers. This prediction is confirmed experimentally in 
Fig. [T] which shows a highly symmetrical G — V curve 
for the n ++ -Si junction. 

If the \E-p- HOMO | energy difference did determine the 
transport barrier (i.e., hole, rather than electron tunnel- 
ing) , then the special undulating J — V behaviour should 
be observed in the n ++ junction for positive applied bias 
on the metal and the p ++ junction should not present 
such behaviour at all. The reason for the latter is that 
for a p ++ -SC - based junction and hole tunneling, 

• with negative bias applied to the metal, holes tun- 
nel from the SC valence band to the metal; 

• with positive bias applied to the metal, holes tunnel 
from the metal to the SC valence band. 

Neither of these flows involves tunneling of holes into 
the SC band gap. Thus, the undulating ln( J) — V be- 
haviour (Fig. 2) and the approximately constant conduc- 
tance (Fig. IT]) for the p ++ -Si-alkyl/Hg junction and lack 
of it in the n ++ one (Fig. [jj) provide direct evidence for 
\E F - LUMO| (or electron^, rather than \E F - HOMO| 
(hole)-dominated transport in both doping types. To put 
this in other words, if transport were HOMO-dominated, 
then the tunneling level and the SC Fermi level would 
be below and above the n ++ -SC gap (at positive bias), 
while both the HOMO and the SC Fermi level would 
be always below the SC gap for p ++ junctions. Thus, 
tunneling via HOMO predicts undulating J — V for n ++ 
junctions, which is not observed experimentally (cf. Fig. 
in the Suppl. Inf. for energy schemes that illustrate these 



alternatives). We note that, for an MIM junction, this 
type of asymmetry cannot be observed, because it origi- 
nates in variation in available density of states for tunnel- 
ing. Indeed, Scott et al have used this exact asymmetry 
and demonstrated that inelastic tunneling is much more 
pronounced in a metal /molecul e-p ++ Si junction than in 
a metal-molecule-metal junction^ 

The conclusion that transport is dominated by the 
\E F — LUMO| energy difference implies that the LUMO 
is closer in energy to the SC Ep for both doping types. 
This can be verified by characterizing the energetics 
of the molecularly modified Si, prior to metal con- 
tact deposition. To this end we use Ultra-violet Pho- 
toemission Spectroscopy (UPS), Inverse Photoemission 
Spectroscopy (IPES), X-ray Photoelectron Spectroscopy 
(XPS) and Contact potential difference (CPD) measure- 
ments. We then study the Si band bending at equilibrium 
after top metal deposition by measuring charge trans- 
port across molecular monolayers (Si-C2H5/Hg) that are 
so thin that they do not present an effective tunnel- 
ing barrier. Finally, we perform a molecular mono- 
layer length-dependent study of SiC„C2n+i/Hg junctions 
where n=14, 16 and 18. These monolayers are thick 
enough to make tunneling the dominant current trans- 
port process. From this length dependence study we ex- 
tract the tunneling barrier height and compare it to the 
barrier that is expected from the UPS/IPES results. 

We begin with determining the band bending (BB) in 
Si in the p ++ - and n ++ -Si-Ci6H33 systems. The position 
of the Fermi level is known in our system from frequently 
repeated UPS and IPES measurements of the Fermi step 
on atomically clean metal (e.g. Au) surfaces. Given the 
nearly 20 A thickness of the SAM layer, the Si valence 
band maximum (VBM) and conduction band minimum 
(CBM) are not visible in UPS and IPES, and we must 
use here the XPS measurement of the position of the 
Si 2p core level, which is clearly detectable through the 
SAM layer, to evaluate the VBM position and the BB. In- 
deed, the energy difference between the core level and the 
VBM, called here (Si2p-VBM) is fixed, and we can use 
the position of either one to determine that of the other. 
Combined UPS/XPS measurements done in our labora- 
tory on carefully prepared hydrogen-terminated Si sur- 
faces give (Si2p-VBM)= 98.9 eV, and we use this value 
here for the determination of the BB. Note that, in spite 
of being obtained with different experimental systems on 
differently prepared surfaces, this value compares well 
with that published in the literatureP^l The Si 2p peaks 
are found at 99.1 eV and 99.8 eV below the Fermi level 
for the p ++ - and n ++ -Si- Ci6H 33 surfaces, respectively. 
Thus, simple considerations give: 

BB p _ si = (E-p — VBM) 

= 99.1 cV- 98.8 eV = 0.2 eV 

BB n _ Si = -Egap — {Ef — VBM) 

= 1.1 cV- (99.8 eV- 98.9 eV) = 0.2 eV 

Both p- ++ and n ++ -Si surfaces are therefore seen to 
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FIG. 3: a) Combined UPS (left of Ef line, blue) and IPES (right to Ef line, red) spectra of Ci6H33-Si(l 1 1) surfaces of 
heavy-doped p ++ - (upper) and n ++ - (lower) Si. The vertical black line marks the experimentally determined Fermi level 
of both samples. The crossings of the tangent lines determine the edges of the monolayer edge-to-edge gap. While the UPS 
measurements yield a single unambiguous edge, which corresponds with the position of the theoretically calculated HOMOp^he 
IPES spectra show two transitions, marking two possible band edges, which are noted as #1 and #2, for both p ++ and n ++ . 
The black vertical arrow marks the position of the estimated tunneling barrier relative to the Fermi level, extracted from length 
dependent J — V measurements of the p ++ Si-alkyl/Hg junction, b) and c) Proposed band diagrams of p ++ - and n++-Si- 
C16H33 surfaces, respectively, which are extracted from the edges of the UPS and IPES results. The Si band lines (solid) are 
slightly bent to reflect the experimentally determined 0.2 eV band bending in both samples and the dashed line marks the 
postion of Fermi level in the Si band gap. The grey shaded areas mark the estimated uncertainties for the positions of the 
actual HOSO and LUSO (see text). 



be slightly depleted with BB= 0.2 eV. The small BB 
on both types of substrate is consistent with the known 
ability of these alkyl monolayers to passivate SiP^I The 
CBM in each case is obtained by adding the 1.1 eV Si 
band gap to the position of the VBM. 

From the position of the photoemission onset of the 
He(I)(21.22 eV) UPS spectra, we find that the work func- 
tions of p++- and n++-Si-Ci 6 H 3 3 are (4.7 ± 0.1) cV and 
(4.2±0.1)eV, respectively. Similarly, the CPD measure- 
ments performed in the dark also yield a 0.5 eV difference 
between the two surfaces, in excellent agreement with the 
photoemission data. Incidentally, CPD also shows these 
values to be independent of the number of carbons in the 
alkyl chain (2,14,16,18) for a given substrate (i.e. p ++ - or 
n ++ -Si).The difference between work functions of p ++ - 
and n ++ -Si-Ci6H 33 (0.5 eV) is smaller than that expected 
from the difference in the Si 2p core level positions for 
the two surfaces (0.7 eV). This may indicates that the Si- 
C bond polarization, which contributes to the interface 
dipole, is affected by the position of the Si Fermi level. 
Strictly speaking, the HOMO and LUMO are concepts 
adapted from a single molecule gas phase model, while 
here we use them for a monolayer in close contact with 
a substrate. Therefore, in the remaining of the text, we 
refer to them as HOSO and LUSO, highest occupied and 
lowest unoccupied system orbital, 40 similar to HOMO* 
and LUMO* in ref. |H 



Fig. [3] presents the combined results of He(II)(40.8 eV) 
UPS and IPES measurmenst both for p ++ (top) and 
n ++ -Si-Ci6H33. We define the edges of the monolayer 
band gap as the crossing between the tangent lines of 
the spectra. While the edge of the UPS spectrum is de- 
fined clearly (both for p ++ and n ++ ), from the IPES 
spectra, two different edges can be defined. The reason 
is that the IPES incident electrons have a larger penetra- 
tion depth than the escape depth of the UPS electrons. 
As a result, the interface states are more pronounced in 
the IPES spectrum than in the UPS spectrum. These 
states, which result from hybridization between the alkyl 
molecular states and the Si VB and CB states at the in- 
terface are known as Induced Density of Interface States 
(IDIS) and were calculated to be localized in the first 
three carbons of the alkyl chain (with respect to the Si 
surface). 38 IDIS is probably the main physical origin for 
the onset of the intensity in the IPES spectrum (Fig. 

As a result, the IPES spectrum has a "two slope" 
behavior, noted as #1 and #2 in Fig. |3]for both the p ++ 
and n ++ samples. We view slope #1 as expressing the 
IDIS and slope #2 originating from the actual LUSO. 

We note that experimentally determined edges of the 
gap do not necessarily correspond exactly to the HOSO 
and LUSO of the monolayer. The onsets of UPS and 
IPES intensity, which define the gap edges, can be af- 
fected by spectral broadening, due to polarization of the 
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molecules in the d ense m onolayer jSHHI by IDIS Poland by 
radiation damage! 44 *^ The HOSO and LUSO positions 
can be found by fitting the calculated density of states 
results, as reported by Segev et al. for the Si-C„H 2n +i 
systenPSJ to the experiment. In their calculations Segev 
et al., identified the actual LUSO to be ~ 0.5 eV above 
the edge, defined by slope #2 . The actual HOSO was 
identified to be ~ 0.3 eV below the UPS edge. The main 
implication of this calculation for the present study is 
that, while there is some experimental uncertainty re- 
garding the actual position of the HOSO, the LUSO and 
band gap width, this uncertainty is unidirectional (i.e., 
the experimental edge to edge gap is the lower limit of 
the actual gap) and a similar situation holds for both 
UPS (HOSO) and IPES (LUSO). 

Figures and c clearly show that for both p ++ and 
n ++ samples, the Fermi level is closer to the onset of the 
IPES spectrum (LUSO) than to the onset of the UPS 
spectrum (HOSO) and the estimated uncertainty of the 
exact position of the HOSO and LUSO is marked by 
a shaded area. The apparent difference in the edge to 
edge gap in Fig(3] for the molecular monolayers on p ++ 
and n ++ , is within this uncertainty. For p ++ -Si-Ci6H33 
samples, the energy difference between the Si Ep and 
the monolayer LUSO (~ 3.2 eV), is smaller than that 
between Ep and the HOSO (~ 4eV). Similarly, for the 
n ++ -Si-Ci6H33 samples, the energy difference between 
the Si E-p and the monolayer LUSO (~ 3eV) is smaller 
than that between the Si Ep and the HOSO (~ 4.4 eV). 

A prominent feature in Fig.[3]is that the intrinsic band 
alignment somewhat varies with the doping density of 
Si. In other words, the energy difference between the 
LUSO (HOSO) and the conduction (valence) band edges 
of the Si is not constant and shifts with the position 
of the Fermi level within the Si band gap. This result 
is surprising, because it is normally assumed that the 
electronic molecular properties (ionization potential and 
electron affinity) are independent of the properties of the 
substrate. Still, this result is highly reproducible experi- 
mentally and the same phenomenon was also observed in 
other hybrid organic-inorganic systems P^l Possibly this is 
related to the above-hypothesised doping-dependent Si-C 
bond polarization. Regardless of this result, which is cur- 
rently under further study, the key experimental finding 
is that, for both p ++ and n ++ Si systems, the Fermi level 
of the system and the carrier transport level in the Si are 
significantly closer to the LUSO than to the HOSO. 

We now ask if the spectroscopic results presented in 
Fig. [3] can efficiently predict the charge transport char- 
acteristics of the alkyl chain monolayer? To answer this, 
we compare the effective tunneling barrier height that 
we can deduce from the J — V measurements with the 
expected barrier that we find from the IPES spectrum 
i.e. \Ep — LUSO| energy difference. To this end we need 
to take under consideration the Si band bending in the 
junction at equilibrium (zero bias), rather than the band 
bending at the Si surface as determined by XPS for the 
free molecularly- modified Si surface, because that band 
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FIG. 4: J — V measurements of p ++ and n++ Si-C 2 H 5 /Hg 
junctions (solid and dashed curves, respectively), compared to 
the J-V behavior of p ++ and n ++ Si-H/Hg junctions (dotted 
curve; p ++ and n ++ curves are identical). The error bars 
represent the standard deviation of multiple measurements 
on 8 different junctions. 

bending can change after forming a metal (Hg) contact on 
the molecules. Because the molecular monolayer is very 
thin (4 — 20 A) and supports a large tunneling current, 
there can be charge transfer between the metal and the 
Si substrate across the molecules, until electronic equi- 
librium is established (equalization of the Fermi levels of 
the Si and the metal). This charge transfer affects the 
Si band bending, but should not change the alignment 
between the Si bands (at the surface) and the monolayer 
HOSO and LUSO, assuming that the molecules do not 
become charged. 

Elsewhere we showed^ that, due do the large surface 
dipole that is introduced by the molecules, depositing Hg 
on a moderately doped n-Si-C n H2 n +i sample results in 
a large band bending in the Si, up to the point where 
the Si surface is in strong inversion and charge trans- 
port in the junction is controlled by minority carrier re- 
combination. For the heavy-doped Si used here, the ef- 
fect of the band bending on transport is usually much 
smaller, because the depletion layer is thinner (~ 10 nm 
compared to ~ 1 \im for moderately doped Si) and al- 
lows some charge tunneling through it. This type of 
charge transport mechanism, which is intermediate be- 
tween thermionic emission and pur e tun neling, is known 
as thermionic field emission (TFE)P2E3 

To isolate the effect of alkyl molecules on the Si band 
bending from tunneling through the molecules, we per- 
formed J — V measurements on a junction with an ul- 
trathin molecular monolayer, C2H5 (Fig. [4J, and com- 
pared the results to those obtained with freshly etched, 
H-terminated Si. Because the C2H5 monolayer is so thin 
(~ 0.5 nm), it presents a negligible charge transport bar- 
rier, compared to the depletion layer in the Si. Combined 
with CPD evidence that the surface dipole is un-altered 
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by molecular length, the ultra-short alkyls allow quanti- 
fying the effect of the depletion layer on net transport. 
The J — V characteristics of the Si-H/Hg junction were 
measured as a control (dotted curve) and found to be 
Ohmic and identical for p ++ and n ++ junctions. The 
results for Si-C^Hs/Hg are very different with the p ++ 
junction exhibiting Ohmic behavior with current density 
even higher than that of the Si-H samples. The n ++ junc- 
tion yields an exponential increase at both bias polarities, 
a behavior that is typical of charge transport via TFEp3 
The high current density of the p ++ junction is consistent 
with our assumption that the tunnel barrier presented by 
the C2H5 monolayer is negligible, and that the dipole di- 
rection works to slightly decrease the Si band bending. 
The same dipole, on the n ++ -Si surface should, therefore, 
increase the band bending. Indeed, when we fit the n ++ 
results at forward bias (positive bias on the metal) to 
an analytical expression for TFE developed by Padovani 
and StrattonpSl we find a barrier height of 0.8 eV. Be- 
cause the Si is heavy-doped and the Fermi level is near 
the CB, this is also roughly the n ++ -Si band bending 
value at equilibrium. 29 Thus, for p ++ -Si junctions there 
is negligible internal barrier in the Si, while such a barrier 
exists for the n ++ -Si junction, where it has a large effect 
on the overall J — V curve. While this does not necessar- 
ily mean that the molecular contribution to transport is 
overwhelmed by band bending in n ++ -Si, as is the case 
with moderately doped Si junctions,^ this band bending 
decreases the current density significantly. Due to the 
large internal barrier in the n ++ -Si, we limit ourselves 
to the molecular length dependent study of the p ++ -Si 
junction, and then only to the data for positive applied 
bias on the metal, where according to Fig. [l] current is 
mostly due to band-to-band tunneling with the density 
of states available for tunneling approximately constant 
with applied bias. The reason is that, as we noted earlier, 
for negative applied bias the dominant charge transport 
mechanism is ISAT, i.e., transport is controlled by the 
recombination rate at the Si surface and cannot be de- 
scribed with conventional tunneling modelsP^ 

Fig. [5] presents the J — V behavior of p ++ -Si-C„H 2n+ i 
/Hg junctions (n=14,16 and 18) (the length dependence 
of the n ++ -Si-C n H 2rl +i /Hg junctions were published 
elsewhere^). It is clear that the qualitative behavior 
of the J-V curve is independent of chain length. It is 
also clearly seen that, as expected for such junctions, the 
current density decreases exponentially with the num- 
ber of methylene units in the alkyl chainP^ This type 
of behaviour is generally described by the Landauer 
relationJ22HH] j n this m odel, th e conductance of a single 
channel (G) is given as! 18 | 52 | 53 l 
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p-HSi-CH 2n+i /Hg (n=14,16,lS) 



G = G c exp {-(31) 



(1) 



where I is the tunnel barrier width, which can be taken 
as the length of the molecule or the thickness of the 
monolayer, /? is the length-decay parameter and G c is 
the contact conductance. From a G vs. I plot at a given 



E 

u 

< 




1.0 -0.5 0.0 0.5 

Metal bias (V) 



1.0 



FIG. 5: Length dependence of current density vs. voltage 
curves for for p ++ Si-C n H 2n +i/Hg junctions (n — 14, 16, and 
18). Data for the Si-H/Hg junctions (dashed) are given as ref- 
erence. The current density exponentially decreases with the 
molecular length. Error bars represent the standard deviation 
of at least 8 different junctions. 



Kipp, we extract both (3 (slope) and G c (intercept). We 
start by discussing the meaning of G c , followed by an in- 
terpretation the extracted /? value. The extracted G c 
value (at V app = 0.1 V) of the p ++ -Si-alkyl/Hg junc- 
tions is 2 • 10~ 4 Go, where Go is the quantum conduc- 
tance (Go = 2q 2 /h = 77.4 ^S). The Gc value of n++- 
Si-alkyl/Hg junctions was 4 • 10~ 7 Gop^ very low com- 
pared to typical Au-alkyl MIM junctions. It reflects the 
large band bending in the Si that was extracted from the 
J —V behaviour of n ++ Si-C2H 5 /Hg observed in Fig. [4] 
The G c found for the p ++ junction is similar to that re- 
ported for an MIM junction with one chemi-contactP^, a 
result that can be rationalized by recalling the near flat 
band conditions in the p ++ -Si-C2H 5 /Hg junction (Fig. 
[4|!^2lTo further verify the correlation between the spec- 
troscopic evidence (Fig. [3| and the transport results, we 
estimate the tunnel barrier (</> t ) from the J — V curves. 
This barrier should, in principle, reflect the energy dif- 
ference between the Fermi level of the electrodes and the 
energy of the nearest allowed state in the interfacial in- 
sulator, as detected by UPS and IPES. In addition, we 
compare the J— V behavior of our p ++ -Si-alkyl/Hg junc- 
tion to previously reported results on the more common 
Au-S-alkyl-S-Au junction. To that effect, we use the gen- 
eral WKB approximation for tunneling that describes the 
relation between f3 and cf>t as 7 : 



(2) 



where /3o is the decay parameter at equilibrium (V^pp — ¥ 
0), h is Planck's constant, q is the elementary charge 
and m e is the electron mass. This is by no means an 
accurate model that reflects the complexity of molecular 
junctions. Nevertheless, as it is the basis for the Sim- 
mons modeP^ that is widely used for the extraction of 4> t 
in molecular junctions, it allows direct comparison with 
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previous results. The difference in dielectric properties 
between vacuum and insulating materials is commonly 
accounted for by using an effective mass for the tunnel- 
ing electron, which is smaller than the free electron mass. 
Using complex band structure DFT calculations, Tom- 
fohr and Sankey calculated the dispersion relation of the 
forbidden band gap of an infinite alkyl chain and showed 
that by using the dispersion relation in vacuum (Eq. [2]) 
with an effective mass of 0.29m G , a lower limit of the 
barrier is extracted.^! Thus, we use this value of effec- 
tive mass to extract an estimated value of <j) t , which we 
can compare to previous studies and to the spectroscopic 
results presented in Fig. [3] 

The average extracted <pt in the low bias range (0 — 

0. 1 V) of the J-V curves in Fig. [H] is (0.84 ± 0.05) A" 1 . 
Inserting this value into Eq. [2| along with a 0.29m c ef- 
fective mass P yields a 2.3 ± 0.3 eV barrier. Having es- 
tablished that \E F — LUSO| dominates charge transport 
(Fig. [lj , we can locate the extracted <j>t value on the IPES 
spectrum (marked by vertical arrow) in Fig. [3] This value 
is in between the IPES edges and, because it is only the 
lower limit of the barrier, we conclude that we have rea- 
sonable agreement with the LUSO of the alkyl monolayer 
(edge #2 in Fig. j3J). A tunneling barrier of ~ 2.3 eV is 
much smaller than the \E F — HOSO| difference of 4eV 
(Fig. [3]) , decreasing further the possibility of significant 
HOSO-contribution to transport. 

The natural question is how far our observation for 
LUSO-dominatcd tunneling in alkyl-Si junctions is ap- 
plicable to alkyl junctions on any substrate? The an- 
swer depends on the position of the molecular levels of 
the alkyl chains relative to the electrode Fermi level. 
UPS/IPES measurements of alkyl thiols on AtPS yield 
a HOSO-LUSO gap of 7.85 eV and significantly smaller 
\E F - LUSO| (~ 3.35 eV) than \E F - HOSO| (~ 4.5 eV) 
(the latter is in good agreement with previou sly repo rted 
\E F -HOSO| for alkyl thiols on Au (~ 5 eV) !^' 56 ' " 1 The 
HOSO-LUSO gap is similar t o the gap measured for alkyl 
chain monolayers on both S i 58 * 59 ! and GaAs^l with dif- 
ferent binding groups (Si-C, Si-O-C and GaAs-P03~C), 

1. e., it is not influenced significantly by binding group and 
substrate type. The range of the edge-to-edge gap in all 
examined samples is ~ (7.1 — 7.8) eV, if the IDIS on- 
set is ignored. These values agree with values reported 
for polyethylene! 5 ^ Thus we conclude that the HOSO- 
LUSO gap of alkyl monolayers, and their position rel- 
ative to the vacuum level is rather independent of the 
contacting electrodes. Therefore, we expect tunneling to 
be dominated by the \E F — LUSO| energy difference for 
any alkyl junctions, except maybe for electrodes of ex- 
tremely high work function (>~ 6eV). 

Our conclusion is at odds with earlier studies on molec- 
ular transport across alkyl chains, which concluded that 
transport is dominated by the \E F — HOSO| barrier, 
especially for the extensiv ely studied Au-(alkythiol or 
alkyldithiol)-Au systemsP ^ 1 23 1 65 ' 66 1 In most cases the 
evidence for hole tunneling is based on transport across 
relatively short alkyl chains, up to C12, and is indirect. 



An exception is the measurement of thermo-power (See- 
beck coefficient), where the positive coefficient was taken 
as direct evidence for hole / HOSO-mediated tunneling.^ 
Remarkably, though, the thermo-power decreases with 
increasing molecular length and, when extrapolated, is 
predicted to change sign (i.e., change from HOSO- to 
LUSO-dominated, hole to electron tunneling) for alkyl 
chains longer than C12. This length-dependence was ra- 
tionalized by invoking Metal Induced Gap States (MIGS) 
as dominating transport, rather than pure molecular 
states. Such an explanation is very reasonable, because 
with the short dithiol molecules used in those experi- 
ments (Au-SC n H 2n+ iS-Au with n = 2, 3, 4, 5, 6, 8) MIGS 
are expected to extend out to several carbon atoms from 
the contacts. 

In early CP-AFM work on short mono and dithiol alkyl 
chains the extrapolated zero molecular length resistance 
increased as the work function of the contacting metal 
decreased and this was taken to indicate an increasing 
tunnel barrier PI Similar to the thermo-power, this effect 
can be due to the MIGS, rather than to alkyl-derived 
energy levels. 

In earlier work with Si as substrate, the type of 
dominant charge carrier in the semiconductor was sug- 
gested to dictate the type of tunneling through (longer) 
alkyl chains. 67 In later work on a system, similar to 
the one studied here, the lack of temperature activa- 
tion was taken as indicating HOSO-mediated /electron 
tunneling^, an effect that may well be present, but is 
overwhelmed by the effects that we have shown here. 
Except for the last two reports, we note that the main 
difference between ours and other experiments is that 
Si-monolayer/Hg junctions (with a several 100 |J.m diam- 
eter contact, rather than a nm-sized one) make it pos- 
sible to measure conductance through longer molecules, 
which allows measuring molecule-, rather than contact- 
dominated transport by reducing the contribution of 
MIGS and IDIS interface effects to net transport. 



III. CONCLUSIONS 

To conclude, the analogy to asymmetry in NDR for 
degenerately-doped MIS of different doping types indi- 
cates that for both p ++ - and n ++ -typc Si-alkyl chain/Hg 
junctions charge transport is controlled by \E F — LUSO| 
(rather than by \E F — HOSO|). Furthermore, by compar- 
ing charge transport through p ++ - and n ++ -Si-C2H5/Hg 
junctions, we showed that there is considerable band 
bending in the n ++ -Si junction, which has a large ef- 
fect on the charge transport (even though the Si is heavy 
doped). The p ++ -Si junction, though, is near flat-band 
in equilibrium, with comparable transport probability to 
that of MIM junctions. Using UPS and IPES measure- 
ments we elucidated the electronic structure of Si-alkyl 
system and showed that \E F - LUSO| < \E F ~ HOSO|, 
in line with the electron tunneling qualitative argu- 
ment. The tunneling barrier height estimated from the 
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WKB approximation also fits well with the IPES-derived 
\E F - LUSO| difference. 
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